SUPERDIFFUSIVITY OF TWO DIMENSIONAL LATTICE GAS 
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Abstract. It was proved Bill.')! that stochastic lattice gas dynamics converge 
to the Navier-Stokes equations in dimension d = 3 in the incompressible limits. 
In particular, the viscosity is finite. We proved that, on the other hand, the 
viscosity for a two dimensional lattice gas model diverges faster than log log t. 
Our argument indicates that the correct divergence rate is (logt) 1 / 2 . This 
problem is closely related to the logarithmic correction of the time decay rate 
for the velocity auto-correlation function of a tagged particle. 



1. Introduction 

It is well-known that although the classical dynamics is time reversible, the 
macroscopic behavior of the fluid, governed by the Navier-Stokes equations, is time 
irreversible. The measure on the time irreversibility is characterized by the viscosity, 
also called the bulk diffusivity of the system. It can be represented as the diffusion 
coefficient of a second class particle in the fluid. Instead of a second class particle, 
one can study the typical behavior of a tagged particle. Once again, even though 
the underlying dynamics is time reversible, the tagged particle is diffusive. The 
diffusion coefficient in this case is called the self diffusion coefficient. The bulk 
and self diffusion coefficients are two different quantities, but they share similar 
qualitative behavior. 

The Green-Kubo formulae represent the bulk or self diffusion coefficients as time 
integrals over the current correlation function or the velocity correlation function. 
In the fundamental work of Alder and Wainwright [Q, it predicts that the time 
decay of the velocity correlation function is of order t~ d / 2 , here d is the dimension 
of the system. Since the decay is 1/t in d = 2 is not integrable, it predicted 
that the self-diffusion coefficient of the two dimensional fluid diverges. In the later 
work, Alder, Wainwright and Gass [5] proposed that the decay in two dimension is 
actually t~ 1 \ogt~ 1 ^ 2 . However, the logarithmic correction cannot be seen from the 
their simulation. This decay rate was also obtained by Forster, Nelson and Stephen 
[S] by a renormalization group method. The later simulation by van der Hoef and 
Frenkel 6 confirmed that there is a discrepancy to the pure algebraic decay t _1 , 
but is far from being able to determine the precise logarithmic correction. 

If we formally integrate the law t~ 1 (\ogt)~ 1 / 2 , we obtain from the Green-Kubo 
formula that the diffusion coefficient diverges as (logi) 1 / 2 in dimension d = 2. For 
higher dimension, the diffusion coefficient is expected to be finite. This was proved 
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rigorously in various settings. For the stochastic lattice gas models considered in 
0], the bulk diffusion coefficient is proved to be finite for d > 3. One key ingredient 
of these lattice gas models is the asymmetric simple exclusion process. For this 
process, both the bulk and self diffusion coefficients were proved to be finite in [T2~] 
and |14j for dimension d > 3. Van Beijeren, Kutner and Spohn [3] predicted via 
the mode-coupling theory that the diffusivity diverges as (logt) 2 / 3 in d — 2 and 
t 1 ' 3 for d = 1. The two dimensional case was proved in |15| : the one dimensional 
case was only partly solved However, related problems in the one dimensional 
case was solved by integrable method [7j. 

In this paper, we shall proved that for the stochastic lattice model, the bulk 
diffusion coefficient diverges in d — 2. Following the method of we derive a 
series of upper and lower bounds to the diffusivity in terms of continuous fractions 
for operators. If we take the first lower bound to the diffusivity, we obtain the 
divergence rate in Theorem l2.f I If one assumes that the dispersion laws of these two 
sequences (upper and lower) of the continuous fractions converge, the divergence law 
(logi) 1 / 2 can be obtained heuristically. See the discussion at the end of this paper. 
Notice that the exponent 1/2 is different from the 2/3 in the case of the asymmetric 
simple exclusion process. Since the lattice gas models are very complicated, we do 
not know if the argument of |15| can be extended to this case. 



We recall the lattice gas models considered in 0] in dimension d. Denote by 
{ej, j = 1, ...,d} the canonical basis of M. d and let £ = {±ei, . . . , ie^}. Let 
V C R d be a finite set representing the possible velocities. On each site of the 
lattice at most one particle for each velocity is allowed. A configuration of particles 
on the lattice is denoted by 77 = {r) x , x £ Z d } where r) x — {rj(x,v), v £ V} and 
ri(x,v) £ {0,1}, x £ Z d , v £ V, is the number of particles with velocity v at x. 

The set of particle configurations is X — ({0, 1} V ) Z ■ 

The dynamics consists of two parts: Asymmetric random walk with exclusion 
among particles of the same velocity and binary collisions between particles of dif- 
ferent velocities. We first describe the random walk part of the dynamics. Particles 
of velocity v perform a continuous time asymmetric random walk with simple exclu- 
sion. A particle at x waits a random, exponentially distributed time then chooses 
a nearby site x + y according to a certain jump law and jumps there as long as the 
site is not occupied by another particle of the same velocity. If there is a particle of 
the same velocity, the jump is suppressed and the particle waits for a new exponen- 
tial time. All particles are doing this simultaneously, and since time is continuous 
ties do not occur. The jump law and waiting time are chosen so that the rate of 
jumping from site x to site x + y is p(y, v) which should be finite range, irreducible 
and have mean velocity v: 



For the sake of concreteness, we take in this paper p(y,v) = unless \y\ — 
£i<i<d \Vi\ = 1 in which case 



for each vector ej and some 7 > large enough for all rates p(e, v) to be non- 
negative. The generator C ex of the random walk part of the dynamics acts on local 



2. The model 




p(±ej,v) = 7 ± (l/2)e,- • v 
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functions / on the configuration space X by 

(£ ex f)(r,) = E^m;*!)!/^)-/!!?)], 

veV xe1 d 
ee£ 

where 

p(x, e, v; rj) — rj{x, v) [1 — r](x + e, v)]p(e, v) 
is the jump rate from x to x + e for particles with velocity v and 

{r](y, v) if w = v and z = x, 
r)(x, v) if w = v and z = y, 
r/(z, w) otherwise. 

The collision part of the dynamics is described as follows. Denote by Q a collision 
set which preserves momentum: 

Q C {(v, w, v', w') e V 4 : v + w = v' + w'} . 

Assume that Q is symmetric in the sense that (v,w, w',v'), (v', w' , v, w), and 
(v',w',w,v) belong to Q as soon as (v, w, v' , w') belongs to Q. Particles of ve- 
locities v and w at the same site collide at rate one and produce two particles of 
velocities v' and w' at that site. The generator C c is therefore 

where the rate p(y,q,r]) is given by 

p(y,q,v) = v(y, v )v(y,w)[i-T]{y,v')][i-i](y,w / )) 

and, for q — (vq, Vi,V2 i Vg,), the configuration r/ y ' q after the collision is defined as 

y q . . _ J T](y, Vj + 2) if z = y and u — Vj for some < j < 3, 
\ z,u ' 1 T](z,u) otherwise, 

where the sum in should be understood modulo 4. 

The generator C of the lattice gas we examine in this article is the superposition 
of the exclusion dynamics with the collisions just introduced: 

c = c ex + c c . 

Let {r?(t) : t > 0} be the Markov process with generator C. 

2.1. The invariant states. We assume that the sets V and Q are chosen in such a 
way that the unique conserved quantities are the local mass Jo and local momentum 
I a , a = 1, . . . ,d: 

loiVx) = ^2v(x,v) , I a (Vx) = ^2(v ■ e a )i](x,v) . 
vev vev 

Examples of sets V and collision dynamics with this property are easy to produce. 
Consider, for example, V = £ and take Q as the set of all vectors (v, w, v' , w') such 
that v + w = v' + w' = 0. An elementary computation shows that the unique 
conserved quantities are total mass and momentum. 0] presents another example 
in d = 3. 
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For each chemical potential A = (Ao, Ai, . . . , Ad), denote by m\ the measure on 
{0, 1} V given by 

1 d 

mx(€) = Z(A) CXP {^ Aa/a(?) } ' (2 - 1} 

where Z(X) is a normalizing constant, 7 (0 = J2 v ev£( v )> MO = X^ e v( u ' e a)£(^) 
for a = 1, . . . , d. Notice that m\ is a product measure on {0, 1} V , i.e., that the 
variables {£(«) : v 6 V} are independent under m\. 

Denote by p\ the product measure on ({0,1} V ) Z with marginals given by 

p x {v ■ r](x,-) = = mx(£) 

for each £ in {0, 1} V and x in Z d . is a product measure in the sense that the 
variables {r](x, v) : x £ Z d ,v 6 V} are independent under p\. A simple computa- 
tion shows that p\ is an invariant state for the Markov process with generator C 
for each A in R + x R d , that the generator C c is symmetric with respect to p\ and 
that C ex has an adjoint £ ex '* in which p(e, v) is replaced by p*(e, v) = p(—e, v). 

The expected value of the density under an invariant state can be computed 
explicitly Fix a vector v in V and define 8 V : R d+1 — > R + as the expected value of 
7](x, v) under py. 

CXp {Ao + Yfa=l X 4 V ■ e a)} 



0„(A) := E^Mw)] = 



1 + exp{A + X) =i ^ a ( v ' 6a )) 



In this formula and below, for a probability measure p, stands for the expecta- 
tion with respect to p. The expectation under the invariant state p\ of the mass 
and momentum are given by 



u„(A) := E^[I a ( Vx )} = Y,{ve a )0vW- 



vev 

(p(A), ui (A), . . . , Ud(X)) is the gradient of the strictly convex function logZ(A). 
In particular, the map which associates the chemical potential A to the vector 
of density and momentum (p,u) — (p, m, . . . , ua) is one to one. Therefore, the 
chemical potential A = (A , . . . , \d) can be expressed in terms of (p, u): there exist 
a subset 21 of M + x R d and functions A a : 21 — > M, a — 0, ... ,d, such that 

A a - A„(p(A),u(A)) (2.2) 

for each A in M + x M. d . Let A = (Ao, . . . , A^). This correspondence permits to 
parameterize the invariant states by the density and the momentum: For each 
(p, u) in 21, let 

v p,u = PA(p.u) ■ 

2.2. Hydrodynamical limit under Euler scaling. In this subsection, we deduce 
the hydrodynamic equation of the system under the assumption of conservation of 
local equilibrium. Fix smooth functions p : R d — > R + , u : R d — > R d . For each 
e > 0, denote by v z p ^ u ^ the product measure on X with marginals given by 

^(.),u(-){ 7 ?( a; '-) =0 = I/ p( e x),u(ex){»?(0,-)=^} 
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for each x in Z d and each configuration £ in {0, 1} V . Assume that particles are 
initially distributed according to V^r.\ u r.y 

For j = 1, . . . , d, let V~ denote the lattice gradient acting on functions / : Z d — > 
K by Vr/(z) — f(z) — f(z — ej). From Ito's formula, for a — 0, ... ,d, we have the 
conservation law 

d 

where Af"(f) are martingales and tojjj are the currents defined by 

d 

CIM = £ ex I a (x,r,) = E V 7<i' 

j'=i 

where 

u %j = J^jlajilx+ej) + /Xeg-v)(ej'v)bx,j(v) and 

b x ,j(v) = r)(x + ej,v)ri(x,v) - (l/2)[r)(x + ej,v) + r)(x,v)] . 

In this computation, the full generator can be replaced by the exclusion operator 
because the collision operator preserves the density and the momentum. Let w° = 
Wo,i for < a < d, 1 < j < d. 

The expectation of the mass current in the j-th direction under the local Gibbs 
state v e , s u ^ is denoted by 7Toj and given by 

n j(p(ex),u(ex)) := E v e^ ^ [w x j] 

= E(«- ej )^(A(p,«)){^(A(p,w))-l} + j(VJ P )(ex) , 
vev 

while the expectation of the momentum currents ir a j are given by 
TT a .j(p(ex),u(ex)) := E„e ()M) [w% tj ] 

= ^2(e a -v)( ej -v)9 v (A(p,u)){e v (A(p,u))-l} + j(Vju a )(ex) . 
vev 

In both formulas, on the right hand side, p and u are evaluated at ex. 

Assuming conservation of local equilibrium, it is not difficult to derive the hy- 
drodynamic equations in the Euler scale for the lattice gas considered in this article 
(cf. [S]). It is given by the system of hyperbolic equations 

dtp + Y,j=i dx^o.j = , 
d t u a + : <), = 0. 

Notice that the factors 7VJ do not survive in the limit due to the presence of a 
second derivative. They will appear, however, in the diffusive scale. 
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2.3. Incompressible limit. Inserting the local equilibrium assumption into the 
conservation laws, on the time scale e~ x t one would obtain the equations 

dtp + ^E^Li^ttoj = 7A/5, 
d t u a + e~ x Yj d j=\ d x^a. 3 = l&u a . 

Consider a system in which the density is near a constant and velocity is of 
order e. Expanding the mass density, the momentum density and the mass and 
momentum currents, we obtain 

u = £U (i) + £ V 2 ) + -- - , 

7 r = 7r< >W 1) +eV 2 > + --- . 

Using these expansions and assuming that the zeroth order terms of 7r are constants, 
we obtain the incompressible Navier-Stokes equations 

d t u { a ] + T,j =1 d Xj ir^ j =^2 itj 1 £ b D^d Xi d Xj uf ) 

for a — 0, . . . , d and uq = p. Here the diffusion coefficient D"'j — jS a ,b^i,j is a 
diagonal matrix. It turns out that this naive computation is correct if the diffusion 
coefficient is instead given by a Green-Kubo formula. 

2.4. Green-Kubo formula. For simplicity, let Aj = for < j < d so that 
0„(A) = 1/2 for every v <G V. Denote this measure by no and let £(x, v) = r](x, v) — 
e v (\) = T 1 {x,v)-l/2. 

Let t x be the shift by x on the lattice, so that (r x r])(y, v) = rj(x + y,v). Denote 
by -c •, • »=-C •, • the scalar product defined on X by 

</,£> = J2 E ^[.f;r x g] 

for two local functions /, g. Here E^ [g;h\ stands for the covariance between g 
and h: E^lg; h] = E^ [gh] — Let Go be the space of local functions 

satisfying 

£ M [5] - and < g, I a > = 

for < a < d. 

Let x — {Xa b, < a, b < d} be the susceptibility which in our context is given 

by 

Xa,b = <C I a , lb > = En [I a (r)o); Ib{r)o)] ■ 
Denote by cr° the part of the current orthogonal to the constants of motion: 

fc=0 

where the coefficients c"' fc are chosen for <r° to belong to C?o- An elementary com- 
putation shows that 
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ex. = (ao, . . . , ad)- Moreover, for 1 < a < d, 

•uSV 

vev 

For simplicity assume that the susceptibility is a constant times the identity: 
Xa.b = K 6 a _b- A straightforward computation shows that this is the case if the set V 
is a cube centered at the origin. Under this assumption, following the computation 
presented in section 2 of ^0], we- obtain that for < a, b < d, 1 < i,j < d, 

= lS a ,b5i,j + -^r- f ds f dr<^at,e rC a^ 

■ttK JQ JO 
^tK Jq Jq 

where 

and, for a local function h, VE Va [h] = {{d/dao)E Va [h\, . . . ,(d/dad)E Va [h]). In 
dimension d > 3, the diffusion coefficients £)"'• (t) converge, as t | oo, to the 
diffusion coefficients given by the incompressible Navier-Stokes equations in 
Subsection 2.3 (cf. 4 

For 9 in R d and r in R d+1 , let 

d d 

D oAt) = E E rJiD^Win ■ 

a,b=0 i,j=l 

We can now state the main result. Let R™ = R n \ {0}. 

Theorem 2.1. Fix in tm<i r in R* +1 . In dimension d — 2, i/iere exists a 
positive constant Cq — Cq{6,t) so that for all sufficiently small A > 0, 

roc 

/ dte- xt tD g , r (t) > CoA-^oglogA" 1 . 

Recall that J °° e~ At f{t)dt ~ A - " as A — » means, in some weak sense, that 
fit) ~ i Q_1 as t — > oo. Theorem l2.1l is therefore stating that the diffusion coefficient 
Dg, r (t) is diverging as log log f in a weak sense. 

Fix a vector in R d , r in R d+1 and let a = ag >r be given by 

d d 

^ = EE r «^°' ( 2 - 3 ) 

a=0j = l 

An elementary computation shows that for every A > 

/>oo 

/ dte~ xt tDg. r (t) = 7A- 2 ||6»|| 2 ||r|| 2 + A~ 2 < ct, (A - C)~ l a > . 
Jo 

Therefore, Theorem 12 . II follows from next estimate on the resolvent. 
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Lemma 2.2. Fix 9 in M.^ and r in There exists a positive constant Co — 

Co(0,r) such that for sufficiently small A > 0, 

•C a, (A - C)- l a > > C loglogA" 1 . 

Notice that the piece of the current 7{/ a (?7ej) — Ia(vo)} vanishes for the inner 
product <C •, • 3>. We may therefore ignore it in the computations below. 

3. Dual representation 

Recall that we are fixing the chemical potentials to be zero: Aj = for < j < d 
and that we denote by /Iq the product invariant measure associated to this chemical 
potential. Unless otherwise stated, < f,g > stands for the inner product of / and 
g in L 2 (jU )- 

Denote by £ the finite subsets of Z d x V and by £ n the subsets of £ with cardinality 
n, for n > 0. For a set A in £ , let ^a be the local function defined by 

^a(v) = n [»/(*,«) -i/2] = n 

By convention "J^ = 1, where </> stands for the empty set. It is not difficult to 
check that {^>a ■ A e £} forms an orthogonal basis of L 2 (/io). In particular, any 
local function / can be written as Xms£ f(^)^A for some finite supported function 
f : £ — ► R. This latter function is called the Fourier coefficients of the local function 
/ and frequently denoted by Tf. A local function / is said to have degree n > if 
(T/)(A) = for any A in £ c n . 

Notice that for any local functions / = X^Aef K^)^A, 9 — Y^Aee 9(^)*a, 

</,5>p„ = £(1/4)" £ f(A)fl(A). 

n>0 Ae£„ 

The factor (1/4) appears because we did not consider an orthonormal basis since 
E tt0 [t(x,v) 2 ] = l/4. 

We say that two finite subsets A, B of Z d x V are equivalent if one is the 
translation of the other. This equivalence relation is denoted by ~ so that A ~ B 
if v4 = B + x for some x in Z d . Let £ n be the quotient of £ n with respect to this 
equivalence relation: £ n = £ n /~, £ = £/~. An elementary computation (cf. [§]) 
gives that 

« f,g »«, = E(V4) n E K^)fl(^) , 

™>! Ae£„ 

where 

1(A) = J2U + z). (3.1) 
Therefore, if we denote by <C •, • ^> the inner product in L 2 (£) defined by 

«f,s»= E(V4r E W)q(a), 

«>! Ae£„ 

where f , g are defined by (|3.1|l , we have that 

=«T/,Tg» . 

The goal of this section is to examine the action of the generators C ex , C c 
on the Fourier coefficients. More precisely, to find operators h ex , L c such that 
L ei T = TC ex , L C T = TC C . 
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3.1. The exclusion operator. We start with the exclusion part of the generator 
which can be decomposed into its symmetric part <S and its antisymmetric part A 
as given by 

V&V x£l, d 3=1 

d 

vev j=i X £Z d 

To examine the action of the symmetric part of the exclusion generator on the 
Fourier coefficients, we first compute S^a ■ An elementary computation shows that 
for each set A in £, 

d 

s^ A = tEEE^.,.-*^ 

veV 3=1 x&L d 

provided A x>x + e , v stands for 

{(A \ {(x,v)}) U {(x + ej,v)} if (x,v) G A and (x + e 3 ,v) £ A, 
(A\{(x + ej,v)})U{(x,v)} if (x + ej,v) e A and (x,v) & A, (3.2) 
A otherwise. 

In particular, for any local function / = J2aKA)^a, a change of variable B = 
A x>x+ejtV gives that 

d 

Sf = 7E*aEE £{f04*,*+e^) - U)} • 

Ae£ v£Vx£Z d j = l 

Therefore, if we define the operator S as 

(sfx.4) = 7 EE p«^,..) - f(^)i , 

we have that TS = ST. 

We turn now to the antisymmetric part. To compute A&a(v)i observe that 
[r){y,w) — 1/2}^ a is equal to ^Au{(y,w)} if (y,w) does not belong to A and is 
equal to (l/4)^ A \{(y,w)} if (y> w ) belongs to A because (r)(y,v) — 1/2) 2 = 1/4. In 
particular, a straightforward computation shows that 

d 

3 = 1 (x,v)eA 
(x+ej ,v)£A 

d 

+ fa •«){(l/4)*A\{(x+e i ,t,)} -*au{M}} • 

j=l (x+ej,v)eA 
(x,v)gA 

In the first term on the right hand side, the second sum is carried over all pairs 
(x, v) in Z d x V, such that (x, v) belongs to A and (x + ej,v) docs not. Therefore, 
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if / = J2ag£ fiA^A is a local function, after elementary changes of variables, we 
obtain that 

Af = £(Af)(A)*A 

provided A = J + + J_ , where 

d 

(J+f)(A) = E E («i'«){f(^\{(» + ei,«)})-f(^\{(a:,«)})} J 
j=i (x,v)eA 

(i+e,,o)£A 
d 

(J_f)(A) = (1/4)^ E fe^KpU^ + e^Hj-pUfet;)})} 

(x+ej ,v)£A 

and T.4 = AT. 

Let If x = § + A. Up to this point we proved that TC ex = II X T. This operator 
h ex , which is not a generator, can thus be decomposed in three pieces, S, J+, J_. 
S is the symmetric part of L e2: and does not change the degree of a function. In 
contrast, J + increases the degree by one, while J_ decreases it by one and J_ is the 
adjoint of — J + : (J_)* = — J+. In particular, J_ + J + is anti-symmetric: 

< Sf,fl > = < f,Sfl > and < J + f, fl > = - < f, J_g > 

so that < (J + + J_)f, g > = - < f, (J+ + JJ_)g » . 

Moreover, a simple computation shows that in L 2 (£) 

Sf = for all functions f of degree one , (3-3) 

J_f = for all functions f of degree two . 

3.2. The collision operator. The remainder of this section is devoted to the 
collision operator. We start defining a generator CI and showing in Lemma 13.11 
below that it is of the same order as C° for our purposes. We conclude the section 
investigating the action of L\ on the Fourier coefficients. 

Fix a site x in Z d and a point q = (v,w, v',w') in the collision set Q. Let 
q' = (v, W, w' , v'). Denote by C x ,q the generator defined by 

(£ x , q f)(v) = p'(x,q,v){lf(v X ' q )-f(ri)} + lf(v X ' q ')-f(v)}} , 

where 

p\x,q,ii) = r)(x,v)r)(x,w)[l-r)(x,v')][l-r)(x,w')] 
+ r)(x, v') T)(x, w') [1 — rj(x, v)} [1 — rj{x, w)] . 
Since the collision set Q is symmetric, we may rewrite the collision generator C c as 

£ c = a/4) E E^- 

zeZ^ q£Q 

Fix a site x in Z d . We start the analysis of the collision operator by examining the 
generator C x ^ q . Since site x is fixed, we omit the index x below so that C q = C Xtq . 
We also denote by C, configurations of {0, 1} V and by £(v) the function — 1/2. 
Since m\ given by (|2.1I) is a product measure on {0,1} V , {IlueB^^) ■ B d V} 
forms an orthogonal basis of L 2 (mo) if IIug^ = ^' wnere stands for the empty 
set. 
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Fix q = (v,w,v' ,w') in Q and let H q = {v, w, v' , w'}. Since C q only changes 
the variables {((u) : u £ H q }, all other variables can be considered as constants so 
that for any subset B of V, 

new = n cm a n cm- 

ueB ueB\H q ueBnH q 

A similar identity holds if we replace £ by £. 
Let 

0i (0 = + eK) - cm - , 

3 (O = £(«KH£K) + £(«)£H£(«>') - £K)£K)£M - £K)£K)£M • 

The index 1 and 3 stand for the degree of the functions involved. Straightforward 
computations give the following identities for degree one functions: 

C q £(v) = £ q £(w) = (1/2)^(0 + 20 3 (e) , 

£ q t(v') = C q t{w') = -(1/2)^(0 - 203(0 • 

Degree two functions vanish under the action of the generator: 

C q Z{ Ul )Z{u 2 ) = 

for ui, U2 in H q . To derive these identities we used that £ q {£(v) + £(i>')} = and 
similar equalities. Degree three functions are such that 

£ 9 £(«)£H£(«') - C q £(v)t(w)£(w') = -(l/8)0i(O - (1/2)03(0 , 
£ q t(v')C(w')C(v) = C q £(v')t(v/)t(w) = (l/8)0i(O + (V2)03(O • 
Finally, degree four functions vanish under the action of the generator: 

A^KMCKKK) = o . 

Here again, to deduce this equality we used that £ q {£{v)£(w){;(v') + £(v')£(w')£(v)} 
vanishes as well as similar identities. 

It follows from the previous formulas that the unique non zero eigenvalue —4 is 
associated to the eigenfunction ip = (pi + 40 3 . In particular, the generator C q can 
be written as 

< ip, Ip > 

where < •, • > stands for the inner product in L 2 (m ). Denote by £ g ,i, C q ,3 the 
operators defined by 

r t , < /, 0i > ^ r , , < /, 03 > , 

£q,lf = -k —L : ^ 01 l £q,3J = -4 03 
< 01,01 > < 03,03 > 

for 03 = 40 3 . Since 0i, 03 are orthogonal, an elementary computation shows that 

-C q < -2£,,i - 2£ 9 , 3 (3.4) 

in the matrix sense. 

If we reintroduce the index x, we obtain the operators 

^ = a/ 4 ) EE £ ^ 
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for j = 1, 3. Notice that both operators keep the degree of local functions. Indeed, 
for a local function /, an elementary computation shows that 

because < 4> x ,q,j, 4>x,q,j >x= 1 for all x, q and j. In this formula < •, • > x stands for 
the inner product with respect to 7710(77(2:, ■)), which means that only the variables 
at site x are integrated and <p x ,q,j are the functions defined by 

<Px,q,i{v) = €(x,v') + £,(x,w') - £{x,v) - £(x,w) , 
^,9,3(7?) = 4£(x, v)£(x, w)£(x, v') + 4£(x,v)£(x,w)£(x,w') 

- 4Z(x,v')£(x,w')Z(x,v) - 4£(x,v')S{x,w')Z(x,w) . 
If / = ^ b for some finite set B, an elementary computation shows that 

< <f>x,q,j >x 4>x,q,j 

has the same degree as ^ b for all x, q and j, which proves the claim. 
It follows from (j3.4j) that 

-L c < -2£l - 2£ c 3 . 

In order to have a tractable algebraic expression for the collision operator, we 
plan to substitute L c by C\ . In order to estimate the third degree terms we use the 
following lemma. 

Lemma 3.1. There exists a finite constant Co such that 

Co' 1 «/,{A-£ ea; -A}-7» 

< «/,{A-£}- 1 /» < Co «/,{A-£ e *-^}-7» 

for every A > and every local function f . 

The proof of this result is based on the lemma below whose proof is similar to 
the one of Lemma 4.2 in [TT] 

Lemma 3.2. Consider a function lo : £3 — » R. Assume that there exists t > 1 such 
that u>((xi, t>i), (x2, V2), (X3, 713)) = if \x\ — x-i \ > £ or \x\ — ^3 > I. Then, there 
exists a finite constant Co depending only on lo such that 

f(x,v)w(x,v) < Co Kx,v)(-Sf)(x,v) 
for every finite supported function f : £3 — > R. 

Proof of Lemma 13.11 We claim that — £§ < — CqS for some finite constant Co- 
Since C% keeps the degree, to prove the claim we only need to show that 

<-cyj>< c Q <-sfj> 

for local functions of a fixed degree. 

Fix 77 > 1 and a local function / of degree n. By definition of £§, taking 
conditional expectations we obtain that 

<-£§/,/> = ]T ]T E[<f,<t> x , q , 3 >l~ . 

q£Qx£Z d 
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Since the velocity set is finite, by definition of 4> x ,q,3, to prove the claim it is enough 
to show that 

>l ] < Co < (-S)f, f > , 

x£l, d 

where B x = {(x,v%), (ar,^), (#, U3)} and V\, v%, V3 are three distinct velocities in 
V. Assume that / = J^B^i-B^B- An elementary computation shows that the 
expectation appearing on the left hand side of the previous inequality is bounded 
above by 

(1/4)" £ f(Bf , 

BZ)B X 

where the sum is performed over all sets B which contain B x . In particular, if for 
a finite set A — {{x\, V\), {x%, V2), (£3, K3)}, we set 

p{Af = (1/4)" J2 « B ) 2 . 

BDA 

where the summation is performed over all sets B which contain A, we just proved 
that 

<-cyj>< J2 E p({(^i),(^2),(x, W3 )}) 2 ■ 

V1,V2,V 3 x£Z d 

By Lemma 13.21 this expression is less than or equal to 

C P{x,v)(-&p)(x,v) < C <-£/,/> • 

Last inequality follows from Schwarz inequality and concludes the proof of the 
claim. Here and below, Co, ao are constants whose value may change from line to 
line. 

Since -C\ < -C S, we have that -C c < -2{C\ + £§) < -C (<S + £?) and 

A - S - C c < a {\ - S - CI} 

for some finite constant ao > 1 and all A > 0. 

On the other hand, since £3 > CqS and since (a — 6) 2 > (1 — e)a 2 — (e _1 — 1)6 2 
for every < e < 1 , a straightforward computation shows that 

-c c > + ^ £ ^^ > -fi^rs + Co^- 1 - 1)5 

for some finite constant Co- Here the factor 1/2 appeared because < ip,ip >= 2 
< 4>i,4>i >= 2 < </>3, ^3 >. If we choose e small enough for Co(e _1 — 1) < 1, it 
follows from this inequality that 

A - S - C c > A - Q-^-Cl - {1-C (e _1 -1)}5 

> Oq^A - S - C\} 

for some finite constant a > 1 and all A > 0. 

Up to this point we proved the existence of a finite constant ao > 1 such that 

a-^A - CI - S} < A - S - C c < a {\ - S - CI} (3.5) 

for all A > 0. 
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It remains to add the asymmetric part of the exclusion generator. Denote by R s 
(resp. R a , R*) the symmetric part (resp. asymmetric part, adjoint) of an operator 
R. It is well known that 

{(R-^y 1 = R*(R S )- 1 R = R s + (R a )*(R s )- 1 R a . 

In particular, for every A > 0, 

({(A-/;)- 1 }") -1 = {X-S-C c ) + A*(X-S -C^A. 

In view of i|3.5|) . there exists a finite constant ao > 1 such that 

ao l U\-S-C\) + A^X-S-CD^A} 

< (X-S-C c ) + A*(X-S - C c )~ l A 

< a {(A-5-£J) + A^X-S-CD^a} 

so that 

av 1 {(X-£ e *-£ c 1 )- i y < {(X-C)- 1 } 3 < a {(X-£ e * -CI)" 1 } 3 , 

which proves the lemma in the case of the inner product of L 2 (hq). The extension 
to the inner product <C •, • is standard (cf. 52)- ^ 
We conclude the section examining the action of £\ on the Fourier coefficients. 
For any local function / = YIaes f(A)^ a, a simple computation shows that 

where 

(Lif)(A) = (l/^^^iA^fiAlUiixyW + fiAlUiixy)}) 

-ftA%U{(x,v)})-ftA%U{{x,w)})} . 

In this formula, q = (v, w, v', «/), A x stands for the set of velocities u such that 
(x,u) belongs to A: A x = {u e V : (x, u) E A} , A x for the set of points (y,u) in 
A with y y^z x: A x — {(y, u) £ A : y ^ x} and if i q (A x ) = 1 if A x is an incoming 
velocity, —1 is A x is an outgoing velocity and otherwise: 

i q (A x ) = 1{A X = {v}} + 1{A X = {w}} - 1{A X = {v'}} - 1{A X = {w'}} . 

With this notation we have that T£f = L^T. Notice that — Lf is a non-negative 
symmetric operator in L 2 (£ n ) and that 

< -Lif, f > = E ^ \ < f> ^-i- 1 >* 



4. Cutoff of large degrees 

For n > 1, let Q n = Ui<k<n£n- Denote by 7r„ the orthogonal projection on 
L 2 (£ n ), by P n the orthogonal projection on L 2 (Q n ) and by C n the operator C ex +C\ 
truncated at level n: C n = P n (C ex + C\)P n . In particular, f = X)n>i 7r "f an( i 
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To investigate the asymptotic behavior of <C cr, (A — C ex — C\)~ 1 cr 3>, for A > 
consider the resolvent equation (A — C ex — C\)u\ = a. In the Fourier space, the 
equation becomes the hierarchy equations 

j;^3UA + (A-§-L^)7r 2 UA = a, 

J^-n-fe+iUA + (A - § - Lf)7r fc u A - JJ+7r fc _iU A = , for fc>3 

because = — J_ and because a has degree 2. The hierarchy starts at degree 
2 instead of 1 because the degree one equation is trivial. Indeed, by Ij3.3|) . (§ + 
Li)7TiUa = 0, J_7t 2 Ua = 0, so that the degree one equation 

-J-7T 2 UA + (A - S - L^TTjUA = 

becomes 7TiUa = 0. Hence 7TiUa plays no role and we can set 7TiUa = 0. 

Notice that we are using the same notation a for the local function defined in 
(|2.3(l and its Fourier transform a : £ — > R which is given by 

d 

<t(A) = 6 j {r (e j -v)+J2r a (e a -v)(e j -v)} (4.1) 

a=l 

if A = {(0, v), (ej, v)} for some » in V, 1 < j < d, and a(A) = 0, otherwise. 
Consider the truncated resolvent equation up to the degree n: 

J+7T3UA,« + (A - § - LJ)7r 2 UA,„ = cr , 

J+7r fe+ iUA : „ + (A - S - L^)7r fc UA : „ - J+7r fc _iUA,„ = 0, 3 < k < n - 1 , 

(A - § - Lf)7T n UA,n ~ J+7Tn-lUA,« = . 

(4.2) 

We can solve the final equation of (|4.2|) by 

7TnliA,n = (A - § - LJ ) ~ 1 JJ + 7T„_ lUA,n • 

Substituting this into the equation of degree n — 1, we have 

7r n -iU A ,„ = [(A - S - Lf) + j;(A - § - L5)- X J + ] J + 7T n _ 2 U A ,„ . 
Solving iteratively we arrive at 

7I"2UA,ri = %i& , 

where the operators {T n ,n > 2} are defined inductively by 

r 2 = (a-s-lj)- 1 , r n+ i = {(a-s-l^ + j;^- 1 ^} -1 . (4.3) 

The truncated equation represents the solution of (A — C n )u\ yn = a and hence 
< ^2^x,nj v >=< cr, (A - C n )~ 1 a > so that 

< cr, (A - £n) -1 £7 > = < Cr, 7^C7 > , 

where, for example, 

|(A-§-L^)+J_(A-§-L^)- 1 J + | , 



3 



T 4 



(A - § - L?) + J_{(A - S - Lf) + J_(A - S - LJ)- 1 ^} J H 
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Lemma 4.1. For each A > 0, <C cr, (A — £2/0+1) x cr 3> is an increasing sequence 
which converges to <C cr, (A — L ex — £l)~ 1 <r 3> and <C cr, (A — C2k)^ 1 <J is a 
decreasing sequence which converges to -C cr, (A — C ex — C\)- 1 <j 3>. 

Proof. Since A — § — is positive, it is easy to show from the definition of the 
sequence of operators T n that < % < T 2 and that T m < T n if T m -\ > T n -\. In 
particular, {72/c, /c > 1} is a decreasing sequence, {T 2 t+i, k > 1} is an increasing 
sequence and 7~2k+i < for any k,j > 1: 

< a, (A - £ 3 rV > < < cr, (A - C^a > < • ■ • (4.4) 
• • • < <C a, (A - £ 4 ) _1 cr > < < cr, (A - £ 2 )~ X cr > . 

To check that <C cr, (A — £ ea: — £^) _1 ct ^ is in fact the limit of these upper and 
lower bounds we use the variational formula. For any matrix M, let M s denote the 
symmetric part (M + M*)/2. The identity {[Af- 1 ^} -1 = M*(M S )- 1 M always 
holds, and thus we have 

< cr, (A - C ex - CI)' 1 a > (4.5) 
= sup {2 < cr, f > - < (A - C ex - £?)f, (A - § - L5)- X (A - £ e * - £?)f > } - 

where the supremum is carried over all finite supported functions f : £ — > R. Note 
that 

« (A - C ex - Cl% (A - § - L5)- X (A - - £?)f » 

= <C f, (A — S — LJ)f > + < Af, (A - S - LJ)- x Af > , 

where A = J + + J_. Hence, 

< cr, (A - £ ea: - LI)- 1 a >= supinf (2 < a - A* , f > 

f fl L 

- <f,(A-S-L5)f > + < fl,(A-S-LJ)fl > } . 

Let a n denote the supremum restricted to finite supported functions in f in L 2 (Q n ), 
and a n denote the infimum restricted to finite supported function q in L 2 (Q n ) so 
that a n T< cr, (A - C ex - CD^a > and a" |< cr, (A - £ e2; - £f)- x ff >. By 
straightforward computation one checks that a n <<S^ cr, (A — £„ + i) _1 cr ^>< a™, 
giving the desired result. □ 

In what follows we will present a general approach to the equations l|4.2|l which, 
from (|4.4|) gives a nontrivial lower bound on the diffusion coefficient. Because it 
gives a sequence of upper and lower bounds, the method has the potential to give 
the full conjectured scaling of the diffusion coefficient. 

5. Removal of hard core 

From Lemma 14. II of the previous section we have a lower bound at degree three. 
However, computations are complicated due to the hard core exclusion. We follow 
[H] to remove the hard core restriction in the formulas and then perform explicit 
computations in Fourier space. By removal of the hard core, we mean replacing 
functions defined on £ n by symmetric functions defined on E n = (Z d x V) n and 
replacing operators acting on E n by operators acting on E n . 
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We first identify a function / : £ n — » R with a symmetric function / : E n — > R. 
Denote by u;„ = (wi, . . . , u„), = (xi, Vi), the points of E n . For n > 1, let 

E n ,i = {w„ : uji ± Wj, for i ± j} (5.1) 

and define 

^ u ^ = |/({wi,...,w„}) if w n G.E„,i, 
1 otherwise . 

With the notation just introduced, 



E Ha I 



[(EfA) 2 ]=^ E «-™) 2 - 



For a function / : E n — > R, we shall use the same symbol (/) to denote the 
expectation 



TZ^ E /K) 



n!4" 

and write the inner product of two functions as (/, g) = (fg) . If / and g vanish on 
the complement of E n ,i, this coincides with the inner product introduced before. 
We also define, as before, <C f,g >= J2 x ez d ( T xf> #)• 

Let E = \J n >iE n , G n = Ui<j< n Ej. We use the same symbol ir n , P n = 
5Zi<j<„ f° r the projection onto E n , G n . As before, there is a simple formula 
for the inner product <C • , • 3>. Consider two finitely supported functions /, 
g : E n — > R. By definition, 

n\ 4™ < /, .g > = ^ /(w n + «)j(w„) = ^ /( w »M w n) , 

zGZ d 

where u„|z= ((xi + z, v\), ...,(x n + z, v n )). 

Denote by ~ the equivalence relation on E n defined by u> n ~ u)' n if for all 
1 < i < n, Vi — v\, Xi — x\ — z for some z in 7L d . Let E n — E n \^. Since summing 
over all sites in E n is the same as summing over all equivalence classes and then 
over all elements of a single class, the previous sum is equal to 

E E + Z ")S(^n + Z) = ^ /( W n)<K W n) 

because f(uj n + z) = /(u> n ). It remains to choose an element of each class. This 
can be done by fixing the last coordinate x n to be zero. In conclusion, 

n!4"«/, ff »=]T /* (wn-i, %>„-i, «), (5.2) 

v£V oj„_iG_E„_i 

where 

/*(w n _i,u) = ^ f{{xi +z,v 1 ),...,(x n _ 1 +z,v n _ 1 ),(z,v)) . (5.3) 

Here again we see that the translations in the inner product effectively reduce the 
degree of a function by one. 



18 



C. LANDIM, J. A. RAMIREZ, H.-T. YAU 



We derive now explicit formulas for the operators §, A+ acting on symmetric 
functions of E n . An elementary computation shows that 

d n 

(s/)(« n ) = tEEE^v^^^./k) 

fc=l i=l i=± 

if u> n belongs to E n< i and (§/)(u> n ) = if uj n does not. Here, for i = ±, 

^k^n = {w 1 ,...,W i -i,(x i + Le k ,Vi),0Ji +1 ,...,0J n ) , 

v^/K) = /(v^w n ) - /( Wn ) . 

Note that S is the discrete Laplacian with Neumann boundary condition on E n> \. 
In the same way, 

n d 

(J+/)(a>„) = ^ ^ y^(efc • t>i)l{:Cj + e fe = gj, « 3 - = «i}V+'/(fa?re) (5-4) 

2—1 J^Z fc— 1 

if u„ belongs to -En,! and (X|_/)(a>„) = otherwise. Here, 

V*+ J 7(a» n ) = /(<)-/(<) 

and the index j in o;^ indicates the absence of uij in the vector oj n : u) 3 n — 
(u>i, . . . , ijJj-i, Wj+i, ■•■,Wn)' Finally, notice that 

n 

(L;/)(w n ) = (1/4) £ £ igfoOlto ^ ^ for fe ^ j} 
<?eQi=i 

[/( (7 .?>"* , n) + /(fj.tu'Wn) — f((Jj :V U} n ) — /((Tj lU ,W n )] , 

if w n belongs to -E n ,i and (Lf /)(<*>„) = otherwise. Here, = (wi, . . . , Wj-1, 

. . .,u n ). 

We now extend the operators S, J+ to symmetric functions not necessarily van- 
ishing on E n> i by formulas analogous to the ones above, except that we drop some 
indicator functions. Let S, J + and LI be the operators defined by: 

n d 

(sf)M = j(af)M = tEEE^W- 

z=l i=± k=l 

n d 

(J+F)(u) n ) = ^2^2^2(e k ■ Vt^ixj + e k = x l ,v i = VjjV 1 ^ F(u; n ) and 

i—1 j^i k—1 

n 

(Llf){u n ) = (1/4) £ £ *,(«,■)!{** # *j for * ^ i} 

96Sj = l 

[/(o-j^/Wn) + f(a JtW 'Lj n ) - f(a JtV uj n ) - f(<Tj, w u) n )} . 

Notice that (LfF) — (J+F) = if < oo and hence the counting measure is 
invariant. Let 

L = S + L C 1 + J+ . 

and denote by L n = P n LP n the restriction of L to G„. Following section 4 in 11 , 
we prove the next result which permits to avoid the hard core interaction of the 
exclusion. 
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Proposition 5.1. In dimension d = 2, there exists a finite constant Co such that 

-^r-Ccr, (A - L n ) _1 o-» < <cr, (A - £„) _1 o-> < C n 4 <cr, (A - L„) _1 cr> . 
for all A > 0. 

The proof of this proposition is similar to the one of Lemma 3.1 in and 
therefore omitted. The main difference is to prove Lemma 4.3 in with S + Lf , 
§ + in place of S, § and this is elementary. 

The special case n = 3 combined with Lemmas 13.11 14.11 tells us that 

< a, (A - C)~ l a > > C < a, (A - L^a > . 

6. Fourier computations 

To bound below <C cr, (A — L^)^ 1 a 3>, define the Fourier transform of a function 
f : (Z d x V)" -> R by 

for p„ £ (R d /27rZ d )™. Here we represented w n — (wi, . ■ ■ , w n ), Wi = (x i: Vi), as 
(x n ,v n ) with x n = (xi, ...,x n )£ (Z d ) n , v n = (vi, . . . ,v n ) £ V™. 

An elementary computation together with (|5.2|) shows that for any local functions 
/, g of degree n, 

</,g> = 77, Un-Ud , An Y\ I f*(Pn-l, V n)S*{Pn-l, V n)dp n _ 1 . 

In this formula, f — Tf, f* is defined by l|5.3|l and f* is the Fourier transform of 
(f*)(x n -i,v n ). Expressing f* in terms of f we further obtain that 

«/,g> = (2n Y n -l)d nHn E /[-*,,]»< T(Pn.»n)fl(Pn,«n)dPn- 

Fix a symmetric function f : (Z d x V)™ — » R. The Fourier transform of the 
discrete Laplacian acting on f is given by 

-Af(p n ,u„) = J(p n ,v n )W(p n ) , 

where 

n n d 

j=l j=l k=l 

if p n = (pi, . . .p n ). Notice that we are using the same notation W(-) for slightly 
different objects. Moreover, for n = 2, a straightforward computation shows that 

d 

J+f(Pl,P2,P3,V 3 ) = -i^^Hv^ = V a2 } (ej ■ V a% ) X 

j=l a 

{sin(ej -p ai ) + sm(e J ■ p„ 2 )}](p<j 1 + p a2 , P* 3 ; v ai , v„ 3 ) , 

where a runs over all permutations of degree three. 

We are now in a position to state the first estimate based on Fourier arguments. 
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Lemma 6.1. Fix a symmetric function f : (Z d x V) 2 — » R. There exists a finite 
constant Co smc/i £/icrf m dimension 2, 

« J+f, (A- A 3 )-V + f » < 
Co [ dpW(p)\log(\ + W( P ))\\%,-p;v 1 ,v 2 )\ 2 . 

The proof of this result is similar to the one of Lemma 3.2 in and therefore 
omitted. 

For < j < 2, let lj = TIj, < j < 2 be the symmetric functions in E n 
associated to the conserved quantities. An elementary computation shows that 

Io(u) = 1, h(v) = ei-v, h{v) = e 2 -v. 

Let Q be the non positive symmetric matrix corresponding to the operator L\ 
acting on functions depending only on one site. Notice that I„, a = 0, 1, 2, are 
eigenvectors of Q. Denote by {I a , 3 < a < |V| — 1} the other eigenvectors of Q 
and by {q a , < a < |V| — 1} the corresponding eigenvalues. Since lo, Ii, I2 are 
associated to conserved quantities, qo = qi = 92 = 0. Let O be the orthogonal 
matrix which diagonalizes Q: 

V = 0*QO . 

For 71 > 1, let 

n 

Q„ = ^I<g>---®Q<8>---<g>I. 

i=l 

Since L\ has an indicator function, we have that < — L c < —Q n . Moreover, Q n 
can be diagonalized by O n = 0® n : V n = 0* n Q n O n , where V n = ® ' ' ' ® v ® 
• • • ® I. Notice that the Laplacian commutes with these matrices. 

As in section0J we can represent (A — Lj,)~ l a in terms of the operators S, J+ 
and L\ to obtain that 

< a, (A - La) -1 * > = <a,{A-A-Lj + j;(A-A - L^)^ 1 J+} 1 a > . 
Since < — L l c < —Q n , the previous scalar product is bounded below by 

< a, {A - A 2 - Q 2 + j;(A - A 3 )- 1 J+} _1 ( t > . 

Recalling that Q 2 = 2 T) 2 2 an d that A commutes with the operators 2 , we 
rewrite the previous expression as 

« o* 2 a, {a - a 2 - v 2 + o;r + (\ - A3)- 1 J + 2 } _1 0*a » . 

To keep notation simple, let = {A - A 2 - T> 2 + 2 J+(A - A3)- 1 J+0 2 }- 1 and 
denote by 7f the projection onto the eigenspace corresponding to the zero eigenvalue 
of T> 2 . By Schwarz inequality, for any function H, 

<ff,Off>> (1/2) < TtH, UttH > - < (1 - W)H, 0(1 - %)H > . 

We claim that -c (1 — it)0 2 <t, 0(1 — tt)0 2 (j 3> is bounded by a finite constant. 
Indeed, by definition of O, 

< (1 - 7f)C*cr, 0(1 - Tf)0* 2 a > < < (1 - ir)0* 2 a, (-D 2 ) _1 (l - n)0* 2 a > . 

Since (1 — ff) is the projection on the positive eigenvalues of T) 2l the previous 
expression is less than or equal to a finite constant which depends on a lower bound 
for the positive eigenvalues of Q and an upper bound for a. 
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The following lemma is needed to estimate <C Tr0 2 a, flnO^cr Recall that 

a = Eo<«<d J2l<j<d ra ®j a j- 

Lemma 6.2. fr0 2 a = if and only if 9 = or r = 0. 

Proof. Assume that Tt0 2 a = and assume without loss of generality that 6\ ^ 
0. Fix the component {0, ei}. Since 7f is the projection of the eigenspace of T> 2 
associated to the zero eigenvalues, n0 2 a — if and only if the scalar product of a 
with (I a , If,) vanishes for < a, b < 2. 

An elementary computation based on the explicit expression 14.1(1 for a shows 
that 

o- •(!„,!(,) = 6 1 ^2l a (v)I b (v){r (e 1 -v)+r 1 (e 1 -v) 2 + r 2 ( ei -v)(e 2 -v)} . 

Setting a = and b = 1 we get that 

0i { r o( e i ' v f + n(ei ■ v) 3 + r 2 (ei ■ v) 2 (e 2 ■ v)} . 
vev 

The last two terms vanish due to the symmetry of V. Hence ro = 0. Repeating the 
same argument with a = b = 1 (resp. a = 1, b = 2), we obtain that r\ = (resp. 
r 2 = 0) . This concludes the proof of the lemma. □ 

We are now in a position to bound below 

< TfO^a, {A - A 2 - V 2 + 0* 2 J* + {\ - A3)- 1 J + 2 }- l TtO* 2 a > . 

Fix a function G in the image of tt. By the variational formula for the 
norm, the previous scalar product with G in place of Tr0 2 a is equal to 

sup{2< G,F > - < F, [A — A 2 - V 2 ]F > - <c C 2 F, j;(A- A3)- 1 J+C 2 F > } , 

where the supremum is performed over all local functions F. Restricting this supre- 
mum to functions in the image of tt we obtain a lower bound. In this case we may 
remove the operator T> 2 because 7f is the projection onto the space associated to 
the zero eigenvalues of T> 2 . 

Lemma lfi.ll and a straightforward computation show that 

< 2 F, j;(a - A3)- 1 J+Q 2 F > 

is bounded above by the right hand side of the statement of Lemma 16.11 Now, 
repeating the arguments presented in the proof of Lemma 3.3 in |llj . we obtain 
that the previous variational formula is bounded below by 

r V f \G(p,-p,V2)\ 2 

°^i h ^ )2 A + 2W(p) + C{W(p)\ log{A + W(p)}| P 

for some finite constants Co, C\. We used here that G is in the image of 7f because 
we needed the function F which maximizes the supremum to be in the image of 7f. 
The factor 2 appeared because W(p, —p) = 2W(p). 

Replace G by n0 2 a. It follows from the previous lemma and elementary compu- 
tations that J2 V \ft®2 cr (Pi ~Pi v 2)\ 2 is bounded below by a positive constant close 
to the origin. The previous expression is therefore bounded below by 



1 

)2 A + 2W( P ) + CiW(p)\ log{A + W(p)}\ 
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for some e > 0. Changing to polar coordinates, since there is no angle dependence, 
we get an integral of the form 

p rdr f p —rdr 



A + r 2 (2-Cilog(A + r 2 )) J Q (A + r 2 ) log(A + r 2 ) 
for some p > 0. By the change of variables u = — log(A + r 2 ) we finally get 

_ r losX du 

J — log p 2 u 

This proves Lemma \2. 21 



log log A 1 



Recall the recursive relation l|4.3|l . If we denote the limit of T n = T. then it 
satisfies the equation 

r = [(A-s-L^ + j^r- 1 ! 

We now assume that the dispersion relation of 7~ _1 is given by 

oo oo 

{E M/ fe)}i i °g( A +{E^)})r 

for some k > 0. Write X^jli W(j3j) = u + W(pi) where u = J2^L 2 W(Pi)- Suppose 
that we are interested in u > A. Then J_T _1 JJ + is approximately given by 

li-e,er U + W(p) + (u + W(p))\\og{u + W(pW ^ ~ l0gU|1 " K 

For u > A we have 

7 1 = (A - § - LS) + J-T" 1 ^ - J-T- 1 ^ 

Thus we have n = 1 — k and this gives the value k = 1/2. This is precisely the 
exponent derived in [21 EI - 

Notice that the exponent is different from the ASEP which takes the value k — 
2/3 ^5]. This is because that the dispersion law is changed only in one direction 
for ASEP. In the fluid model considered here, the collision operator spread the 
dispersion law to all direction. This causes the exponent on the right side of (|f> . 1|) 
to be 1 — K, compared with 1 — n/2 in JS]. This sketch is certainly not rigorous 
since it is not even clear where the operator T should be defined. Though one can 
try to prove upper and lower bound as in |15|. it is not clear that the off-diagonal 
terms can be controlled. However, the exponent 1/2 seems to be very convincing. 
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